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ABSTRACT 

We present Herschel observations of the core of the Perseus cluster of galaxies. Especially in- 
triguing is the network of filaments that surround the brightest cluster galaxy, NGC 1275, 
previously imaged extensively in Ha and CO. In this work, we report detections of far- 
infrared (FIR) lines, in particular, [Cn] 158 fim, [Oi] 63 //m, [Nn] 122 /mi, [Oib] 145 /im 
and [Om] 88 fim, with Herschel. All lines are spatially extended, except [Om], with the [Cn] 
line emission extending up to 25 kpc from the core. [Cn] emission is found to be cospatial 
with Ha and CO. Furthermore, [Cn] shows a similar velocity distribution to CO, which has 
been shown in previous studies to display a close association with the Ha kinematics. The 
spatial and kinematical correlation among [Cn], Ha and CO gives us confidence to model the 
different components of the gas with a common heating model. 

With the help of FIR continuum Herschel measurements, together with a suite of coeval 
radio, submm and infrared data from other observatories, we performed a spectral energy 
distribution fitting of NGC 1275 using a model that contains contributions from dust emission 
as well as synchrotron AGN emission. This has allowed us to accurately estimate the dust 
parameters. The data indicate a low dust emissivity index, ft « 1, a total dust mass close to 
10 7 M G , a cold dust component with temperature 38 + 2 K and a warm dust component with 
temperature of 1 16 + 9 K. The FIR-derived star formation rate (SFR) is 24 ± 1 M Q yr _1 , which 
is in agreement with the FUV-derived SFR in the core, determined after applying corrections 
for both Galactic and internal reddening. The total infrared luminosity in the range 8 /im to 
1000/imis inferred to be 1.5xl0 n L Q , making NGC 1275 aluminous infrared galaxy (LIRG). 

We investigated in detail the source of the Herschel FIR and Ha emissions emerging from 
a core region 4 kpc in radius. Based on simulations conducted using the radiative transfer code, 
cloudy, a heating model comprising old and young stellar populations is sufficient to explain 
these observations. The optical line ratios indicate that there may be a need for a second 
heating component. However, stellar photoionization seems to be the dominant mechanism. 

We have also detected [Cn] in three well-studied regions of the filaments. Herschel, with 
its superior sensitivity to FIR emission, can detect far colder atomic gas than previous stud- 
ies. We find a [Oi]/[Cii] ratio about 1 dex smaller than predicted by the otherwise functional 
Ferland (2009) model. That study considered optically thin emission from a small cell of gas 
and by design did not consider the effects of reasonable column densities. The line ratio sug- 
gests that the lines are optically thick, as is typical of galactic PDRs, and implies that there 
is a large reservoir of cold atomic gas. This was not included in previous inventories of the 
filament mass and may represent a significant component. 
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1 INTRODUCTION 

The Perseus cluster of galaxies is preeminent among the class of 
cool-core galaxy clusters (those with gas cooling times shorter than 
the Hubble time). This is in large part due to its close proximity 
(z = 0.01756), allowing detailed studies to be conducted in varying 
astrophysical contexts. It is the X-ray brightest galaxy cluster and 
has a strongly peaked surface-brightness profile. The intracluster 
gas in the inner few tens of kpc has a very short radiative cool- 
ing time (200 Myr to 300 Myr). In the absence of heating, the ex- 
pected X-ray mass deposition rate is several 100 M Q yr _1 . FUSE 
observations, on the other hand, suggest an actual cooling rate of 
~ 30 Mq yr 1 jBregman et al.|2006| > and XMM-Newton RGS obser- 
vations suggest an even lower residual cooling rate of 20 M Q yr~' 
( |Fabian et al.|2006| and references therein). 

Perseus is the prototype of cluster radio "mini-haloes" (e.g. 
|Pedlar et al.|1990| |Gitti et al.|2002) . The brightest cluster galaxy 
of Perseus (a giant cD galaxy), NGC 1275, is host to a powerful 
radio source, 3C 84. It has a Seyfert-like spectrum and a bolo 
metric radio core-luminosity of order 10 43 erg s~ ! (e.g|Vermeulen 



|et al.| |"l994). 3C 84 is inferred to undergo episodic bursts of ac- 
tivity, blowing out jets of plasma which interact with the intra- 
cluster medium. Chandra X-ray images, when compared with ra- 
dio emission at different wavelengths (e.g. Bohringer et al.|1993| 
Fabian et al.|[2000l [20061 1 Sanders et al.||2005HSanders & Fabian| 



2007| > reveal the extent to which the central radio source has caused 



havoc in the intracluster-medium (ICM) in the form of bubbles of 
various kinds (inner, outer, ghost), sound waves and shocks. On 
VLBI/VLBA (milliarcsec) scales, the radio morphology consists of 
a bright southern jet (e.g. Vermeulen et al. 1994, Taylor et al. 2006) 
believed to be pointed towards the observer and a dim northern 
counter jet pointed away. The radio source shows a clear double- 
lobed morphology on smaller scales (~ 80") that extends into a 
more amorphous one on larger scales (~ 170"). 

One of the most intriguing aspects of Perseus is the spec- 
tacular network of ionized (Ha) and molecular (CO and Ht) gas 
filaments well beyond the optical stellar emission of the bright- 
est cluster galaxy (BCG). A significant number of BCGs in cool- 
core (CC) clusters show similarly extended filamentary structures, 
such as NGC 4696 (Centaurus), Abell 1795 and Hydra-A jjohn" 



stone et al.|1987l|Heckman et al.|1989||Sparks et al.|1989|[Craw- 



ford et al. 1999, McDonald et al. 2010). The source of excitation of 
the filaments is currently one of the most pertinent issues in our un- 
derstanding of CC galaxy clusters. In the case of NGC 1275, the fil- 
aments extend as far as ~ 50 kpc out from the cluster-centric AGN, 
with the mean surface-brightness declining much slower than the 
inverse-square law. Photoionization from a central AGN thus seems 
unlikely. Similarly, ionizing radiation from hot stars, such as O and 
B type, has been ruled out based on anomalous emission lines in the 
spectra of the filaments ( Johnstone & Fabian 1988 , Johnstone et al. 
|2007|[Ferland et a l. 2008). Motivated by the observations of strong 
molecular hydrogen lines in NGC 4696 and NGC 1275, |Ferland] 
|et al.| (2009) showed that non-radiative heating, such as collisional 
heating from ionizing particles, can produce the observed emis- 
sion. The importance of collisional excitation by energetic (ioniz- 
ing) particles was suggested more than two decades ago by |John-| 
|stone & Fabian| (1988). Candidate sources for these particles are 
either cosmic rays or the thermal electrons of the X-ray emitting 
intracluster medium. Motivated by the spatial correspondence be- 
tween the brightest low-energy X-rays and the Ha filaments, Fabian 
|et al.| |20TTJ considered the penetration of cold filaments by the 
surrounding hot X-ray gas through reconnection diffusion. More 



recently, |Sparks et al.| (2012[ > reported a detection with the HST 
ACS camera and also the COS spectrograph, of [Ov]/U549 line 
emission spatially coincident with the Ha line emission in M87 in 
Virgo. The [Civ]/11549 line emission is indicative of gas at ~ 10 5 K. 
They suggest the origin of this line emission as being due to ther- 
mal conduction, i.e., the transport of energy by hot electrons from 
the hot ICM to the cold filament gas. 

Various independent studies have provided strong evidence of 
the presence of dust in cool-core BCGs. These include dust con- 
tinuum observatio ns (e.g.|Edge et al.|1999||Chapman et al.|2002| 
|Egami et al.|2006|[OTJea et al.|2008[|Rawle et al.|2012| > and HST 
observations of BCGs with dust absorption features in them (e.g. 
McNa mara et al.|1996||Pinkney et al.|T996l|Laine et al.|2003||Oonk| 
|et al.||2011| l. Furthermore, observations of H 2 and CO molecular 
gas (e.g. |Donahue et aLpOOOj pdgel[200T] |Edge et al.|[2002l |Sa^ 
|lome et al.|2006[ > suggest that there are substantial amounts of dust 
present which provide shielding. In M ittal et aTjpOl l[ >, we demon- 
strated the need for an overabundance of dust (low gas-to-dust mass 
ratio) and metallicity to explain measurements in NGC 4696. In 
most of the cases, the gas-to-dust mass ratios are consistent with 
Galactic values ( Sparks et al. 1989, Edge 2001 1, however, there is 
clearly an overabundance of dust in the filaments relative to the sur- 
rounding hot gas. The origin of the dust in the filaments is not yet 
clear. More recently, Voit & Donahue (2011} have suggested that 
the main source of dust in BCGs of CC clusters may be the stars 
of the central galaxy themselves. Thus, a stellar origin is plausi- 
ble although in the case of NGC 1275 that seems unlikely given 
the filaments exist out to a large radius from the galaxy core. Even 
though the sputtering timescale for dust grains is short compared 
to the gas cooling time, optical observations have revealed mag- 
nificent dust lanes in some objects, which in many cases correlate 



spatially with the Ho- filaments (Crawford et al.|2 005 ; Sparks et al. 
fl^89l|Don"ahue & Voit| 1993} . |Fabian| ( | 1994| > and |Voit & Donahue 
(19951 have proposed that gas-phase reactions create dust in the 
cold clouds that cool out of the ICM, which can then lead to an 
increase in dust through accretion onto existing grains. In contrast, 
several studies argue (e.g. |Sparks et al.|1989| |Farage et al.|2010| > 
that the ionized gas and dust filaments originate from the stripping 
of a dust-rich neighbouring galaxy which may be in the process of 
merging with the BCG. The origin of the molecular and ionized 
filaments mixed with dust is to date an open question. 

In this paper, we present Herschel observations of the core of 
the Perseus cluster. The main goal of this work is to investigate the 
source of the various emissions originating from the filaments. Far- 
infrared data from Herschel ( Pilbratt e t al.|2010| l have proven to be 
useful diagnostics of the heating mechanisms that account for the 
filaments in cool-core clusters (e.g. Rawle et al. 2012 Mitt al et aT| 
[20TT]|Edge et al.|2010a|b||Pereira et al.|2010) . This work is part of a 
Herschel Open Time Key Project (PI: Edge) aimed at understand- 
ing the origin of cold gas and dust in a representative sample of 
1 1 BCGs. We describe the data used and the analysis in section[2] 
We present some basic results in section [3] and move on to dis- 
cuss the heating mechanisms prevailing in the core of NGC 1275 
in section [4] and those prevailing in the filaments of NGC 1275 
in section [5] We finally give our conclusions in section [6] We as- 
sume throughout this paper the ACDM concordance Universe, with 
H = 71 hn km s _1 Mpc" 1 , £l ra = 0.27 and fi A = 0.73. This 
translates into a physical scale of 1" = 0.352 kpc and a luminosity 
distance of 75.3 Mpc at the redshift of NGC 1275 (z = 0.01756). 
This distance is consistent with the independent distance inferred 
from the 2005 Type la supernova SN2005mz ( |Hicken et al.|2009) . 
The right ascension and declination coordinates in figures are in 
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J2000 equinox. Lastly, NGC 1275 comprises two systems - a low- 
velocity system (LVS) consisting of gas at 5200 km s -1 associated 
with the BCG and a high-velocity system (HVS) consisting of gas 
at 8200 km s -1 associated with a foreground galaxy north-west of 
the BCG. In this work, NGC 1275 refers to the LVS system only 
unless otherwise mentioned. For calculation of line velocities we 
assume the velocity of the LVS to be the systemic velocity. 



ScanMap observing mode, particularly designed to detect extended 
emission (> 100"). The pipeline employs a second-level deglitch- 
ing algorithm, which uses the redundancy in a pixel to flag outliers, 
so that bright sources are not erroneously flagged as glitches. For 
SPIRE photometry, the data were reduced using the pipeline for the 
LargeScanMap observing mode and the naive map-maker. 



2.2 Hey data 



2 DATA AND ANALYSIS 



2.1 Herschel data 



We used the PACS spectrometer (Po glitsch et al.|2010) to observe 
the [Cn] line at 157.74 pm and the [Oi] line at 63.18 fim, the two 
primary coolants of the interstellar-medium (ISM). The [Cn] and 
[Oi] fine-structure lines are very often the brightest emission lines 
in galaxy spectra. In addition, we observed [Oib] at 145.52 yum, [Sii] 
at 68.470 fim, [Nn] at 121.90 yimi and [Ora] at 88.36 fjm. The [Oi] 
and [Cn] lines were observed in the raster-mapping mode, consist- 
ing of 5 raster lines and 5 points per line with a line step of 23.5", 
whereas the rest of the lines were observed in a single-pointing 
mode. The observational parameters are summarized in Table [T] 
The line observations were conducted in the line-spectroscopy 
mode using the chopping and nodding technique (using a chopper 
throw of 6') to subtract the telescope background, the sky back- 
ground and the dark current. 

The PACS photometric observations were made in large-scan 
mapping mode at a speed of 20"s~' at blue-short (BS) (70 yimi), 
blue-long (BL) (100 /an) and red (R) (160 /mi) wavelengths (PI: 
E. Sturm, ObsIDs: 1342204217, 1342204218, 1342216022 and 
1342216023). The scans consisted of 18 scan line legs of 4' length 
and of a cross-scan step of 15". The "scan" and orthogonal "cross- 
scan" observations were individually calibrated before being com- 
bined into a single map of 9' x 9'. The PACS photometer has a 
resolution of 5.2", 7.7" and 12" at 70 /mi, 100 /mi and 160 /mi re- 
spectively. The PACS photometer performs observations at BS and 
BL simultaneously with the R band so we have two sets of scans 
in the R band. The SPIRE photometric observations was also made 
in the large-scan mapping mode and the data were recorded simul- 
taneously at 250 /mi, 350 /mi and 500 /mi (PI: E. Sturm, ObsID: 
1342203614). The SPIRE photometer ( |Griffin et al.1[2010| has a 
resolution of about 18", 25" and 36" at these wavelengths, respec- 
tively. 

The basic calibration of the data (spectral and photomet- 
ric) was done using the Herschel Interactive Processing Environ- 
ment (HIPE) jOtt|2010| > version 7.0 CIB 1931. For the PACS spec- 
tral data, the standard pipeline routines described in the PACS data 
reduction guideline (PDRG) were adopted to process the spectral 
data from their raw to a fully-calibrated level. HIPE 7.0.1931 con- 
tains PACS calibration files that provide the response calibration 
based on in-orbit measurements. Hence, no ground-to-flight cor- 
rection factors had to be applied. The PACS cubes were rebinned in 
wavelength using the Nyquist-Shannon sampling, corresponding to 
oversample=2 and upsample=l. The spatial full-width at half max- 
imum (FWHM) varies from 8" for the [Sii] 68.47 /mi line to 11" 
for the [Cn] 157.74 /mi line. The line fluxes were determined us- 
ing the method described in Mittal et al. (201 1). Briefly, the routine 
specproject was used to obtain a final projection of the different 
pointings and nods onto the sky plane. These maps can be read- 
ily used to conduct 'aperture photometry' and measure fluxes. For 
PACS photometry, the data were reduced using the pipeline for the 



The Ho- flux measurements were made using the continuum- 
subtracted data from |Conselice et aT (200 1). Based on comparison 
with our newer HST data ( F abian et al.|2008| >, the calibration of the 
WIYN data appears to over-estimate the flux by a factor of about 
three. The source of this discrepancy is yet not clear and is under 
thorough investigation (Johnstone et al. 2012, in prep.). To verify 
the calibration offset of the WIYN continuum-subtracted image, we 
calculated the total flux in the WIYN Ha image in counts s -1 and 



converted it into erg s cm - using the conversion given in Con- 



|selice et al.|<f200T>. We were only able to derive the total luminosity 
quoted in |Conselice et al.| (|2001), and which is in good agreement 
with |Heckman et al.| ( |1989} , after lowering the total flux by a fac- 
tor of three. Hence we have scaled down the measured WIYN Ha 
fluxes by a factor of three. 

The HST broad band filter F625W admits light from the 
[Oi]i6300, [Nn]A6583 and [Sn]A6731 doublets as well as Ha. The 
ratios of these lines to Ho- are variable with position in the nebula. 
Such spectroscopic data are available for only a small fraction of 
the total area covered by the entire nebula. For the WIYN data the 
contamination is only from the [Nii]/16583 doublet. This is the main 
reason why we prefer to use the WIYN data. 

|Hatch et al.| {2006[ > found a radial gradient in the ratio 
[Nn]/16583/Ha, which may be due to a spatially varying metallic- 
ity or excitation mechanism (see also Johnstone & Fabian 1988| >. 
They used the multi-object spectrograph instrument on Gemini and 
found that the [Nn]/16583/Ha ratio varies from 0.5 to 0.85 in the 
Horseshoe region. Assuming that the [Nii]/16583/Ha ratio does not 
change significantly with azimuthal angle, we subtracted the con- 
tribution of [Nn]/16583 from the measured Ha flux in the Horse- 
shoe, southwest and Blue Loop knots, using an average ratio of 
[Nn]/i6583/Ha=0.65. The Gemini measurement close to the nu- 
cleus indicates a [Nn]/16583/Ha ratio close to unity, therefore we 
halved the measured Ha flux in the core region. Similarly, we also 
corrected for the [Nn]/16548 line, usually a third in intensity of the 
[Nn]^6583 line Patch et al.|2006) . 

Due to the uncertainties in the [Nn]/16583/Ha ratio and the 
calibration of the WIYN data, we caution the reader concerning the 
absolute values of the Ha fluxes. However, the relative values of Ha 
flux between spatial positions should be accurate. The Ha ampli- 
tude plays a role in section]?] where the Herschel and Ha measure- 
ments are used as constraints to determine the heating mechanisms 
giving rise to the various emissions. Fortunately, there are enough 
FIR-derived constraints that the best-fit model parameters do not 
rely solely on the Ha flux measurement. 



2.3 Dust Extinction 

SCUBA observations of NGC 1275 have been used to infer the 



presence of a large amount of dust (6x 10 7 M Q ) present i Irwin et al. 
2001). The source of dust is not yet clear and understanding it is one 
of the goals of this study. However, dust also presents a hindrance 
due to the extinction it causes at optical and higher frequencies. 
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Table 1. Herschel PACS spectroscopy observational log of NGC 1275 at a redshift of 0.01756. All the lines were observed in the line spectroscopy mode and 
on the same day: 30 lh Dec. 2009. 



Line 


Peak Rest A 


ObsID 


Duration 


Bandwidth 


Spectral FWHM 


Spatial FWHM 


Mode 




(/an) 




(s) 


(jim) 


(kms^ 1 ) 


(/«n) 


(kms- 1 ) 






OI 


63.184 


1342189962 


9600 


0.266 


1250 


0.017 


79 


9.4" 


5x5 raster, step size 23.5" 


CII 


157.741 


1342214362 


8600 


1.499 


2820 


0.126 


237 


11.1" 


5x5 raster, step size 23.5" 


NO 


121.90 


1342214363 


3440 


1.717 


4180 


0.116 


280 


10.6" 


pointed 


Olb 


145.525 


1342202581 


3440 


1.576 


3215 


0.123 


250 


9.7" 


pointed 


OIII 


88.356 


1342214363 


3680 


0.495 


1660 


0.033 


110 


8.5" 


pointed 


Sil 


68.473 


1342202581 


3840 


0.218 


945 


0.014 


62 


8.3" 


pointed 



Velocity [km/s] (CII line) 
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Wavelength [urn] 



Velocity [km/s] (OI 63 //mline) 
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Velocity [km/s] (Olb 145 //m line) 
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.4.0 64.1 64.2 64.3 64.4 64.5 64 

Wavelength {/im 



147.4 147.6 147.8 148.0 148.2 148.4 148.6 148.8 

Wavelength [z/m] 



Velocity [km/s] (Nil line) 

-1500 -1000 -500 500 



A 



1000 1500 




123.6 123.8 124.0 124.2 124.4 124.6 

Wavelength [am] 



Velocity [km/s] (OIII line) 

-1500 -1000 -500 500 1000 1500 




89.4 89.5 89 



9.8 89.9 90.0 90.1 90.2 90.3 90.4 

Wavelength ]//m] 



Figure 1. The forbidden far-infrared line detections in the center of NGC 1275 made with the Herschel PACS instrument. The lines are spatially integrated. 
Top Row: [Cn] 157.74 pm (left), [Oi] 63.18 pm (middle) and [Om] 145.52 pm (right). Bottom Row: [Nil] 121.90 /im (left) and [Om] 88.36 pm (right). 



In this work, both the Galactic and internal extinction corrections, 
such as for Ho- (section |2T2] l and far-ultraviolet measurements (sec- 
tion |4.3fr , were calculated with the help of the mean extinction laws 
given in Cardelli et al. ( 1989 ). We assumed the Galactic extinction 
law,R v = 3.1, and nnE(B-V) value ofO.163 from the NASA/IPAC 
Extragalactic Database (NEEQ. Internal extinction was calculated 
using the observed Balmer decrements, which were compared to 
the Case-B value of Ha[HB= 2.86. For the filaments located far out 
from the core, an internal extinction of E(B — V) = 0.38 was es- 
timated based on the Galactic-extinction corrected Balmer decre- 
ment, Ha[HB= 4.2, as measured in the Horseshoe knot (Ferland 
|et al.|[2009} . For the core, a similar internal extinction was esti- 
mated, E(B- V) = 0.37, based on the Galactic-extinction corrected 
Balmer decrement, Ha-/Hy8=4.08, as measured 18" SW of the nu- 
cleus by |Kent & Sargent| ( [T979l >. Note that |Kent & Sargent|p979) 
obtain an internal reddening of E(B- V) = 0.43 assuming a Galactic 



http://nedwww.ipac.caltech.edu 



extinction of E(B - V) = 0. 1. We obtain a lower internal reddening 
due to the higher Galactic extinction adopted. 

The correction factor for the internal reddening is model- 
dependent. If the Ho- and HB emissions are produced by another 
mechanism than case-B recombination, then the Balmer decrement 
will be different from the expected value of 2.86. For example, if 
particle heating is responsible for the Ha and HB emissions, then 
the intrinsic Balmer decrement will be higher and consequently the 
deduced internal reddening lower. 



3 RESULTS 

Of the six lines observed, we detected all except [Sii] at 68.47 pm. 
Even though only [Cn] and [Oi] observations were designed to de- 
tect extended emission, all the detected lines except [Om] are spa- 
tially extended. The integrated line profiles are shown in Fig.[T]and 
their relative spatial extensions in Fig. [2] where the pixel thresh- 
old has been set to SNR>2. The SNR corresponds to the ratio of 
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Table 2. Integrated line properties. Also given is the 3-<x upper-limit for the [Sii] line flux. The spatial extents are based on visual inspection. 



Line 


A dim) 


Offset (k 


ms" 1 ) 


FWHM (km s" 1 ) 


Line Flux 


Spatial Extent (Radius) 






z bcg 




Obs. 


Intrinsic 


(l(r 18 W/m 2 ) 


(") 


(kpc) 


[Oi] 


64.298±0.002 


39±9 


-61±9 


383±20 


375±20 


2525.8±73.7 


34" 


12 


[Cn] 


160.510±0.002 


-0.6±4 


-100±4 


419±9 


347+1 1 


2205.3±26.5 


71" 


25 


[Nit] 


124.03 1±0.004 


-24±10 


-124+10 


558±22 


482±25 


125.0±2.8 


28" 


10 


[OlB] 


148.091±0.003 


32±7 


-68±7 


458±17 


384±20 


150.5+3.1 


29" 


10 


[Om] 


89.900±0.005 


-37±17 


137+17 


375±42 


358±44 


65.6±4.1 


< 9.4" 


< 3.3 
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the line peak to the standard deviation of the data about the fitted 
model. Listed in Table|2]are the integrated line properties. 

The [Cn] line is detected over a spatial region of ex- 
tent 140" (50 kpc). Shown in the left panel of Fig. [3] is a 
Ha+[Nn]/i6583, 6548 image displaying the ionized gas filaments 
taken with the WIYN 3.5 m telescope ( |Conselice et al.||200i) . 
Shown in the right panel of Fig.[3]is the FIR [Cn] emission with the 
pixel detection threshold set to SNR > 1. The middle panel displays 
the Ha emission smoothed to match the resolution of the [Cn] line 
using the ciao tool 'aconvolve' and the pyraf tool 'blkavg', with the 
[Cn] line contours overlaid (in red). The [Cn] emission traces the 
Ho- emission very well despite a much lower resolution. Both [Cn] 
and Ha reveal a central elongation, about 20 kpc in total extent, 
with an east-west alignment. This elongation is clearly visible also 
in the recent narrow-band imaging of a ro-vibrational transition 
line of molecular hydrogen (H 2 ) ( |Lim et aL||2012) , which shows 
very good overall morphological resemblance with Ho- emission. 
In |Mittal et aklpOll) , we showed a close spatial and kinematical 
correspondence between different emissions, such as Ha and [Cn], 
in NGC 4696, the BCG of the Centaurus cluster. NGC 1275 shows 
a similarly tight correlation, both spatial and kinematical, between 
[Cn], Ha and CO (section|3TTJ, suggesting a common heating pro- 
cess of the gas. 

There are three well-studied regions in the Ho- filaments, seen 
also in the [Cn] map, which we describe briefly below: (1) The 
Horseshoe is a filament to the north-west of the center of BCG 
( |Conselice et al.|20 01 1. The filament of gas appears to originate in 
the core of the galaxy, rising a projected distance of 26 kpc through 
the ICM before looping around towards the center again. (2) The 
southwest knots, at a projected distance of ~ 23 kpc from the cen- 
ter of the BCG, are a part of the southern filaments. 3) The knots 
in the southeast filament, also known as the 'Blue Loop', are at a 
projected distance of ~ 16 kpc. The loop is so called due to the 
blue optical colour of the knots this region, indicating these are star 
forming sites. Canning et al. (20101 studied this region in detail and 
found a star formation rate of ~ 20 M Q yr . We refer to the knots 
as the 'Blue Loop' knots from now on. All three regions have been 
studied in detail previously at several wavelengths ( Conselice et al.| 



2001 



2011 



Salome et al. 2006 2008, Canning et al. 2010, Salome et al. 
Hatch et al.|2005||2006||Fabian et al.|2003||2008||Johnstone| 



etal. |20 07, Feiiand et al. 2009|>. 

The three regions are marked in Fig. [5] We extracted the [Cn] 
line spectra from them and the core region (see below). The spec- 
tra are shown in Fig. [4] The Horseshoe region is weaker in [Cn] 
emission than the SW and Blue Loop knots. This is consistent with 
the study of Sal ome et al.| ( [20iT[ l, where they find this region has 
weaker CO emission in comparison with other filaments. There is 
no detection of [Oi] in any of the three regions of the filaments. 
The line strengths and ratios in these regions and the core provide 
strong constraints for studying the excitation mechanism(s) in dif- 
ferent parts of the BCG. This is followed-up in detail in section|5] 

Brauher et al.| ( |2008| > report a flux of (2.5 ± 0.2) x 10 ' 5 W nr 2 
for [Oi] and (1.2 ± 0.1) x 10" 15 W nT 2 for [Cn], based on obser- 
vations with the Long Wavelength Spectrometer (LWS) on the In- 
frared Space Observatory (ISO). While the [Oi] flux measurement 
of (2.53 ± 0.07) x 10~ 15 W nT 2 obtained in this study compares 
well with the ISO measurement, our [Cn] flux of (2.21 ± 0.03) x 
10~ 15 W nr 2 is higher by a factor of two. The reason for this dis- 
crepancy is that Brauher et al. (2008 1 obtain the [Cn] flux based on 
the assumption that NGC 1275 is a point source. The FWHM of 
the LWS is 75". However, [Cn] emission clearly extends beyond 
the FWHM of the LWS (~ 140" in diameter). For this reason, we 



believe that Brauher et al. (2008) underestimate the line flux in [Cn] 
by about a factor of two. As a rough check of this hypothesis, we 
convolved our [Cn] map with a gaussian with a FWHM of 75" and 
obtained a flux of 0.97 x 10~ 15 W irr 2 within a 75" aperture diame- 
ter. Despite the fact that this test uses a simple gaussian rather than 
the true ISO LWS beam profile, the reduction of measured flux is 
consistent with the result of Bra uher et al.| (2008 ). 

3.1 Kinematics 

Fig. [5] the [Cn] velocity distribution shows the full extent of the 
complex kinematical structure of the filaments in NGC 1275. The 
colour scheme is such that the red shaded regions represent red- 
shifted gas (with positive velocities with respect to the systemic 
velocity of the BCG), whereas green and blue shaded regions rep- 
resent blueshifted gas (with negative velocities with respect to the 
systemic velocity of the BCG). The velocity pattern is likely to be 
the combination of inflowing and outflowing gas, along with pro- 
jection and small-scale rotation effects (see below). 

The dashed circle at the center marks the 'core' region (see 
Table 7}, wh ich shows a gradient in the [Cn] line velocity. |Wilman| 
|et al. " (2005) conducted near-infrared spectroscopy of the cen- 
tral region in NGC 1275 with the UIST IFU instrument on the 
United Kingdom Infrared Telescope. They detected ro-vibrational 
H 2 emission originating from a region ~ 50 pc from the nucleus, 
with the results indicating a strong velocity gradient in the peak 
position of the H 2 line. The measured rotation curve in the central 
IFU slit placed east-west along the nucleus shows a velocity change 
from 150 km s east of the nucleus to -100 km s west of the nu- 
cleus. From the sharp decrease in the magnitude of the velocity on 
either side of the nucleus, they concluded that the molecular gas 
is distributed in a disk-like structure with the rotation axis oriented 
north-south. This is also consistent with the double-horn structure 
seen in the CO(2-l) spectrum extracted from the same core region 
(Salome et al. 201 1 1, which may be interpreted as indirect evidence 
for a central rotating disk. The resolution of the [Cn] data is not 
good enough to resolve any rotational structure over the scales ob- 
served by Wilman et al. ( 2005 1. The smallest structure in [Cn] that 
can be resolved is about 4 kpc. Similarly, the spectral resolution of 
the [Cn] data (230 km s~') is much poorer than that of the CO(2- 
1) line observed with the IRAM 30 m telescope on Pico Veleta 
(40 km s~'). From Fig. [5] the rotational pattern claimed by | Wilman 



|et al.| ( (2005] > and[Salome et al. (20fT|l is not surprising. However, it 
is difficult to say whether the double horn feature is due to small- 
scale disk-like rotation or large-scale flows. If there is a rotating 
disk present in the center, it involves a relatively small gas disk of 
~ 5 kpc radius. The majority of the gas distribution sampled by 
[Cn] does not show any regular rotation. 

Also interesting is the redshifted ridge of gas passing through 
the center of the BCG with a north-south extension. The CO veloc- 
ity measurements made by Salo me et al.| < |20 1 1| > in regions 4 and 
21 marked in Fig. 1 of their paper provide a confirmation of the 
presence of this redshifted gas (also see Lim et al. 2008 ). This wide 
vertical distribution of gas has a cylindrical symmetry and is sug- 
gestive of material being dragged upward and downward. 

The upper left panel of Fig.|5]shows that there is a negative ve- 
locity region on either side of the center along the east-west direc- 
tion. The east region extends all the way along to the south, where 
it merges with the Blue Loop knots. The blueshifted components 
on both sides of the major axis of the emission are also visible in 
the kinematics inferred from the CO spectra (Sa lome et al.| 2006l. 
The CO flux and velocity distribution shows a close association 
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Figure 3. : Left The continuum-subtracted WIYN Ha line emission with a spatial resolution of 0.6" to 0.8" jConselice et al.| 2001 1. Middle: The Ha image 
smoothed to the resolution of [On], with the [Cn] line contours overlaid in red. Right: The Herschel [Cn] line emission with a spatial resolution of 11" (this 
work). The optical Ho- emission and infrared [Cn] emission show a remarkably similar distribution. Marked are the Horseshoe knot in the northwest 20 kpc 
from the BCG, the SW1/SW2 knots in the southwest 21 kpc from the BCG and the Blue Loop knots in the southeast 16 kpc from the BCG. 



with the optical filaments seen in Ho- (e.g. Salo me et al.|2006| >. In 
Fig. [6] we show the kinematics of the gas in the inner 90"x90" 
in three phases - optical, FIR and millimeter - as represented by 
Ha [taken from |Conselice et al. 1(2001) ], [Cn] and CO(2-l) [taken 
from |Salome et a l. ( 2006)] line emissions. This plot shows a clear 
overlap in the redshifted gas, with absolute velocities greater than 
5264 km s -1 (red symbols), along the ridge associated with all three 
line emissions. Similarly, there is an overlap in the blueshifted gas, 
with absolute velocities less than 5264 km s (blue symbols), west 
of the ridge. Especially remarkable is the curvature in the ridge 
seen in both Ha (red filled squares) and [Cn] (red crosses). From 
this we conclude a kinematical correlation among [Cn], Ha and CO 
emissions, which further reinforces the idea that different emissions 
have the same origin. 

Radio observations of central radio sources often provide use- 
ful insights in understanding the kinematics of cool-core BCGs. 
Overlaid on the bottom panels of Fig. [5] are contours of the ra- 
dio emission associated with 3C 84 at 1.4 GHz and 74 MHz. The 
bottom left panel shows 1.4 GHz radio emission at two different 
resolutions - 2" (white contours; courtesy of VLA/NRACQ and 
5" (yellow; kindly provided by Greg Taylor). The higher resolu- 
tion contours reflect a north-south radio morphology, demonstrat- 
ing a good alignment with the redshifted [Cn] contours. The lower 
resolution contours reflect an inverted-S shaped morphology, pre- 
viously noted in several studies (e.g. Pedlar et al. 1990, Bohringer 
|et al.|19 93 , Fabian et al. 2000). Although the southern jet shows a 
correlation with the [Cn] emission, the tip of the northern jet does 
not. The 74 MHz contours (also kindly provided by Greg Taylor), 
shown in yellow in the bottom right panel, indicate a reversal of 
the east-west component of the jet direction on both sides of the 
core, wherein the northern jet once again overlaps with the [Cn] 
emission. Conselice et al. ( 2001) noted a similar alignment between 
the linear extension of the Ha filaments to the north and the low- 
brightness radio emission seen at low frequencies. 

Several studies have presented scenarios wherein radio out- 
bursts are responsible for the dredge-up of cold, metal-rich gas 



2 The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 



from the core in the direction of the buoyantly rising radio plasma 
(e.g. |Simionescu et al.|[2008| [2010) |Gitti et al.|20TT] |Revaz et al.| 



|2008[ 1 Tremblay et al.||2012a|b| >. Evidence for these scenarios is 
based on the X-ray-derived temperature and metallicity maps, 
which show a spatial correlation between radio emission and cool 
gas extending away from the core with a metal content higher than 
that of the ambient medium. The positive correlation between the 
radio emission and the redshifted ridge of [Cn] gas in NGC 1275 is 
reminiscent of cold gas being dredged up by the radio lobes. On the 
other hand, this interpretation implies that both the radio jets (jet 
and counter-jet) are pointed away from the line-of-sight. On mil- 
liarcsec scale, the radio morphology comprises a one-sided jet, such 
that the southern jet is deemed to be approaching, and the north- 
ern jet (counter-jet) receding ( e.g. |Pedlar et al. 1990| |Vermeulen| 
|et al.|1994[rDhawan et al.|1998[|Taylor et al.|2006| >. Both the radio 
jets reveal complex kinematics in the plane of the sky. The south- 
ern jet is initially elongated along position angle (PA) ~ 235° (at 
~ 10 mas) but suddenly bends toward PA~ 160° (at ~ 20 mas) and 
continues though a series of such bends, at approximately the same 
angles, out to about an arcminute (Pedlar et al. 1990, Dhawan et al. 
1998 ). The northern jet exhibits a similar complex structure. On the 
kiloparsec scale, it is therefore possible that, similar to the kinks 
observed in the plane of the sky, the jets undergo bends along the 
line-of-sight, such that both the jets are receding. There are dredge- 
up interpretations offered for NGC 1275, suggesting that Ha gas 
is being dredged up by the radio source. |Fabian et ah]p003") rec- 
ognized two filaments in the north-west, including the Horseshoe, 
bent on either side of the north-west ghost bubble. They showed 
that the Ha emission associated with the Horseshoe is just behind 
the bubble and is likely dragged out by it. Similarly, Sanders et al. 
( 2005 1 discovered a high-abundance ridge using Chandra observa- 
tions, which they hypothesized is formed by material entrained by 
a fossil radio bubble. 

Lastly, the [Cn] velocity structure may be related to the dis- 
ordered motion of gas clouds at larger radii that is not affected by 
the inner radio structure. Shown in Fig.|7jis the [Cn] linewidth in 
km s . An interesting feature of this map is that one of the re- 
gions the linewidth peaks is adjacent to the redshifted ridge where 
the line velocity flips from positive to negative on the eastern side 
of the galaxy. This is where the CO(3-2) HARP maps shows dou- 
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Figure 4. The [Cn] emission from the Horseshoe knot in the northwest (upper left panel), the SW1/SW2 knots in the southwest (upper right panel) and the 
Blue Loop knots in the southeast (lower left panel). These apertures are marked in Fig. [5] We also show the [Cn] spectrum for the core region (lower right 
panel), marked in the upper left panel of Fig. [5] 



ble lines implying that there are multiple gas components along the 
line of sight (Edge et al„ in prep.). However, the linewidths in this 
region may also be high due to the superposition of rapidly varying 
line-of-sight velocity elements. 



3.2 Spatial variation of [Oi] to [Cn] ratio 

Here we briefly discuss the spatial variation of [Oi] to [Cn] ratio and 
what may be inferred from it. The resolution of [Cn] line is slightly 
poorer than the [Oi] line. For compatibility, we convolved the [Oi] 
line map with a gaussian such that the resulting FWHM is similar 
to that of the [Cn] line. Fig.[8]is obtained by dividing the resulting 
smeared [Oi] map by the [Cn] map (the threshold for both [Cn] 
and [Oi] was set to SNR>2). This figure indicates that while [Oi] 
is stronger than [Cn] in the core, the opposite is true at larger radii, 
namely, the [Cn] line emission becomes stronger at radii larger than 
~4 kpc, which corresponds to the spatial resolution of [Cn]. 

The relative strength of [Oi] to [Cn] is an indicator of gas 
density, such that a higher [Oi]/[Cn] ratio represents a higher den- 
sity gas for the reasons described in section |4~T) Hence a higher 
[Oi]/[Cn] ratio in the center implies a relatively denser gas in the 
cluster core, as expected. 



3.3 Dust SED 

We detected emission at all three PACS wavelengths and all three 
SPIRE wavelengths. The PACS images are shown in Fig. [9] Dust 
emission can usually be modelled as a simple modified black body 
function; NGC 1275, though, poses a complication. This is because 
of the strong radio source at the center of NGC 1275, with a large 
contribution in the sub-millimeter and FIR range (Irwin etal.|2"00T) . 
The radio source, additionally, shows a large-amplitude variabil- 
ity on timescales of decades (see Nagai et al. 2012, and references 
therein). Detailed monitoring at 3 mm has shown an increase in flux 
from 3.5 Jy in 2002 to 1 1 Jy in 2010 < |Trippe et al.pfJTT) . 

In order to obtain accurate dust parameters, it is essential to es- 
timate the fraction of the FIR flux originating from the synchrotron 
emission from 3C 84. In the following, we attempt to simultane- 
ously fit the dust and AGN emissions. We list in Table [3] and show 
in Fig. [10] the various flux measurements (black filled diamonds) 
that were used to obtain an optimal model. While the dust contri- 
bution to the SED is constrained by the PACS (blue open diamond) 
and SPIRE (red filled triangles) data, the synchrotron contribution 
is constrained by radio and sub-mm data. The two sub-sets of sub- 
mm flux-densities clustered around 1 mm and 3 mm (yellow filled 
squares) correspond to the IRAM Plateau de Bure Interferometer 
(Trippe et al. 2011) observations. These measurements are from a 
monitoring program and chosen so as to be closest in time (span- 
ning days between 2010-08-21 to 2010-08-29) to the SPIRE and 
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Figure 5. Upper Left: The detailed velocity distribution of the [Cn] gas. The colour-bar represents the velocities in km s . The dashed circle represents gas 
compatible with rotation. Upper Right: The velocity contours superimposed on the Ha image. The dashed blue contours correspond to blueshifted (negative) 
velocities and the solid red contours correspond to the redshifted (positive) velocities. Bottom Left: 1.4 GHz radio contours at two different resolutions, 
2" (white) and 5" yellow, overlaid on the Ha image. Bottom Right: 74 MHz radio contours at 25" resolution (yellow) overlaid on the Ha image. 



PACS observations. The 5 GHz, 8 GHz and 15 GHz radio mea- 
surements (purple diagonal-plus symbols) correspond to the obser- 
vations made in 2010 with the 26-m telescope at the University 
of Michigan Radio Astronomy Observatory (UMRAO, courtesy of 
M. Aller and H. Aller). We used only the 15 GHz data point for the 
fitting since the spectrum turns over below this frequency, indica- 
tive of synchrotron self-absorption. The 50 /im and 25 fim measure- 
ments (green filled circles) correspond to the IRAS observations 
(Mosh lFet al.|1990) . To obtain a robust estimate of the tempera- 
ture and mass of the plausible second (warm) dust component, we 
also used the 20 yum and 30 fim continuum data (red open crossed- 
squares) determined from the Spitzer IRS observations ( |Weedman| 
|et al.||20"05] >. We did not fit 6 /urn, 10 fim and 15 fim data due to 
a possibly increasing contribution from a passively evolving pop- 



ulation of stars at wavelengths < 15 fim. The AGN synchrotron 
contribution to the total emission is expected to be small at wave- 
lengths shorter than 70 /jm and so the variation in the flux density 
due to the different times of the IRAS/Spitzer observations from 
the rest of the fitted data can be neglected. 

To model both the components - the dust emission and AGN 
synchrotron emission - we used the following fitting function: 



S v — S syn,F + S dust,v > where 



S o(^) ' if V < Vbrcak 



^dust,v — 



(l+z) 3 



[B r (T d ) - B r (T cmb )] (l-e-^ M «>), 



(1) 

(2) 

(3) 
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Figure 6. The kinematics of the three phases, optical, far-infrared and mil- 
limeter, represented by Ha, [Cn] and CO(2-l) line emissions. Show are 
the redshifted (red symbols) and blueshifted (blue symbols) velocities for 
the three phases. The [Cn] velocities are shown as plus (blueshifted) and 
cross (redshifted) symbols, the Ho- velocities are shown as filled circles 
(blueshifted) and squares (redshifted) and CO(2-l) velocities are shown as 
open down-triangles (blueshifted) and up-triangles (redshifted). 




Figure 7. The [Cn] linewidth in km s~ 1 . The dashed blue contours corre- 
spond to blueshifted (negative) [Cn] velocities and the solid red contours 
correspond to the redshifted (positive) [Cn] velocities (same as the upper 
right panel of Fig. [5}. 



Eqn.[2]quantifies the synchrotron emission assumed to have a bro- 
ken powerlaw form. So denotes the normalization at y =100 GHz 
(3 mm), a t and a 2 and the (negative) powerlaw indices on either 
side of the powerlaw break frequency, Vb rea k, respectively (an rep- 
resenting the powerlaw at radio frequencies and cr 2 representing the 
powerlaw at submm/FIR frequencies). We find that for NGC 1275 
it is not necessary to modify the synchrotron emission with an ex- 
ponential term like |Privon et al.|(2012| > did to model the spectrum of 
Cygnus A. The need for including an exponential term in the case 



Table 3. A compilation of the fitted radio, sub-mm and IR flux-densities 
for NGC 1275. The columns are (1) wavelength (2) instrument (3) year 
of the observation (4) aperture (available only for this work) and (5) the 
measured flux-density. Note that the measurements correspond to the total 
flux-densities for the given instrument. 



A (fim) 


Instrument 


Year 


Aperture 


Flux (mJy) 


20 


SPITZER IRS d 


2004 


- 


2410±241 


25 


IRAS C 


1983 


- 


3539±176 


30 


SPITZER IRS d 


2004 




382()±382 


60 


IRAS C 


1983 




7146±286 


70 


PACS Herschel a 


2011 


55" 


7405±741 


100 


PACS Herschef 


2010 


55" 


8541±854 


160 


PACS Herschel 3 


2010 


55" 


6979±1396 


250 


SPIRE Herschel" 


2010 


60" 


3805±571 


350 


SPIRE Herschel" 


2010 


70" 


3095±464 


500 


SPIRE Herschel" 


2010 


112" 


2992±449 


1153 


IRAM PdB b 


2010 




7220±1083 


1303 


IRAM PdB b 


2010 




7741±1161 


1428 


IRAM PdB b 


2010 




7299±1095 


2710 


IRAM PdB b 


2010 




10000±1500 


2913 


IRAM PdB b 


2010 




10713+1607 


3019 


IRAM PdB b 


2010 




10330±1550 


20675 


UMRAO e 


2010 




23230±130 



a This work. The errorbars correspond to the absolute flux uncertain- 
ties: 10% at PACS BS and BL and 20% at R, and 15% at all SPIRE 
wavelengths; b |Trippe et ak]j201 1) ; c |Moshir et"aL]jl990t ; d |Weedman] 
et al. 1 2005 1; e University of Michigan Radio Astronomy Observatory 
data for year 2010 (courtesy of M. Aller and H. Aller). 



of Cygnus-A arises from the lack of synchrotron emission beyond 
the submm regime. NGC 1275, on the other hand, shows evidence 
of an AGN contribution all the way out to submm and, possibly, 
IR wavelengths ( Krabbe et al. 2000 1. Hence, we preferred a broken 
powerlaw, such that \a 2 \ > to an exponential term to represent 
the slowly decaying synchrotron emission. 

Assuming the dust to be in thermal equilibrium, Eqn. [^quan- 
tifies the dust emission, which is a blackbody function modified by 
a term that depends on the dust optical depth, r v , defined as 



Din' 



(4) 



The modification from a standard blackbody function is due to the 
fact that for typical dust temperatures, the dust grains are smaller 
than the peak wavelength of the Planck function and hence do not 
radiate as perfectly as a blackbody. Here, M d is the dust mass, r d 
is the dust temperature, k v = K vo {vlvaf is the dust absorption coef- 
ficient and we adopted k vq = 1 m 2 kg -1 at 1200 GHz (250 fim) 
(Hildebrand 1983). f} is the dust emissivity index, which based 
on empirical results likely lies in the range 1 to 2. B V (T) is the 
Planck function at frequency v, and temperature T. B v (T cmb ) is the 
contribution from the cosmic microwave background at T cmb = 
2.73(1 + z) K. CI is the solid angle subtended by the source, here 
assumed to be the total extent of the FIR emission at 70 /urn. 

The powerlaw indices for the AGN emission, a { and a 2 , were 
constrained to have negative values, so that the synchrotron emis- 
sion decreases with increasing frequency. Dunne & Eales ( 2001 ) 
showed that a two-component dust model with B close to 2 better 
fits the observed SEDs than a single-dust component. This is also 
consistent with our findings based on the SEDs of cool-core BCGs 
investigated so far ( |Edge et aLl|2010a|b| |Mittal et aTT||201 1^ . We 
therefore fitted the data with a model comprising a cold and warm 
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Table 4. The best-fitting model parameters of the function describing the FIR emission (Eqn[TJ for the dust emissivity index, ft = 1 . For comparison, we also 
give the parameters for ft = 1 .5 and /3 = 2. The uncertainties correspond to the 1-cr errorbars. 



Parameter 


/?= 1.0 


P= 1.5 


P = 2.0 


Cold dust temperature, c 

Cold dust mass, c 
Warm dust temperature, T± w 

Warm dust mass, Mj w 
Synch, powerlaw norm., So 
Synch, powerlaw index 1, al 
Synch, powerlaw index 2, al 
Synch, powerlaw break frequency, i-Wak 


(38.0±2.0) K 
I0 6.96±0.12 Mo 

(115.5±8.5)K 

1() 4.71 ± 0.17 Mo 

(10434±393)mJy 
-0.41±0.02 
-1.08±0.28 

(2.6±0.5)xl0 n Hz 


(35.0±1.6) K 

10 6.91±0.12 Mo 

(107.0±7.4) K 

10 4.37 ± 0.17 Mo 

(10364±393) mjy 
-0.42±0.02 
-0.92±0.24 

(2.6±0.6)xlO n Hz 


(31.8±1.3) K 

10 6.89 ± 0.12 Mo 

(95.7±6.0) K 

10 4.60 ± 0.18 Mo 

(10448±393) mjy 
-0.41±0.02 
-0.83±0.21 

(2.5±0.7)xlO u Hz 


Derived Quantities 


Dust luminosity (8 /an- 1000 /mi), Lfir 
Star formation rate, SFR 
Total gas-to-dust mass ratio 


(1.4±0.05)xl() u L Q 
24±1 M Q yr 1 
4500 to 7800 


(1.3±0.05)X10 U Lq 
23±1 M a yr" 1 
5000 to 8000 
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22±1 M Q yr -1 
5300 to 9200 
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Figure 8. [Oi] to [Cn] ratio. [Oi] has been convolved with a gaussian such 
that the resulting full-width at half maximum (FWHM) equals that of the 
PSF of the PACS spectrometer at the wavelength of [Cn] . The circle shown 
in the bottom left corner represents the size of the FWHM. The ratio map 
shows that [Oi] is stronger than [Cn] in the core. The trend reverses at ~ 
4 kpc from the center. 



dust component parametrized by temperature and mass (T i c , M d c ) 
and (7d_ w , Md, w ), respectively. We used the Levenberg-Marquardt 
non-linear least square fitting algorithm from Numerical Recipes 
to obtain the best-fit model parameters appearing in Eqn. [T] M dc , 
Td.c, M& iW , T i w , /3, So, cm, ct2 and v br eak- We explored a range of 
values between 0.5 and 2.5 for the dust emissivity index,/?. The^ 2 - 
minimization gave a best-fit value /? which was < 1 . Note that there 
is a strong degeneracy between /3 and the dust temperature such 
that all values of /3 explored yield models, with varying tempera- 
tures, that are compatible with the observed SED. I.e., for any of 
the explored values of ft between 0.5 and 2.5, the best-fit \ 2 value 
was less than 11; for a chi-squared distribution with 17 - 8 = 9 
degrees of freedom (because we are not optimizing over /?), this 
corresponds to a p-value (Gregory 2005) greater than 0.10, i.e., a 
model consistent with the data. Based on the work of IDunne &l 
|Eales| ( (2"001[ >, however, p is expected to lie between 1 and 2, and 



so we fixed /3 to unity and found the best-fit values for the other 
parameters subject to that choice. 

The best-fit model is shown in Fig. [To] (black solid line), as 
are the AGN (thin-dashed line) and the dust contributions (thick- 
dashed lines; red representing the cold dust component and blue 
representing the warm dust component). In addition to the fitted 
flux-densities, we show the supplementary Spitzer/IRAS and UM- 
RAO data, along with SCUBA Jlrwin et al.||200T] > and WMAP 
( |Wright et al.|2009) data from around year 2000. The SCUBA (big 
magenta crosses) and WMAP data (small indigo crosses) clearly 
demonstrate the strong variability in the radio source and, hence, 
the need for coeval flux-densities for SED fitting. Similarly, the 
Planck (gray plus symbols Planck Collaboration et al.|20TT| > and 
UMRAO data (open green squares) from the year 2009 fall slightly 
below the IRAM and UMRAO data, both from the year 2010, re- 
spectively. 

The best-fit parameter values were used to calculate the to- 
tal dust emissivity by integrating the emission between 8 /(m and 
1000 /im. Assuming that young (< 10 s yrs), hot stars dominate 
the interstellar radiation field across the UV-optical band, the star 
formation rate can be estimated using the Kennicutt relation ( |Ken-| 
|nicutt|[T998[ >. In early-type galaxies, including BCGs, the cooler 
emissions (> 100 /mi) may arise from dust heated by a passively 
evolving old stellar population, which warrants caution in the Ken- 
nicutt calibration. The best-fit parameters for/3 = 1 and the derived 
quantities are given in Table|4](we also give the best-fit parameters 
for p" = 1 .5 and /3 = 2.0 for comparison) and the predicted dust and 
AGN flux contributions are given in Table [5] A total gas-to-dust 
mass ratio between 4500 and 7800 was estimated using the molec- 
4 x 10 10 M Q derived in[i 



Salome et al. 



(2006). 



The 



ular gas mass 

total gas mass depends upon the conversion factor used to calculate 
the atomic plus molecular mass. We used a factor of 1.36 as given 
in Edge (2001). The total gas-to-duss mass ratio, although high, is 
within the range of the derived mass ratios in other cool-core BCGs 
( |Edge|2001^ . At the other extreme is NGC 4696, for which a 3-cr 
upper-limit of ~ 450 was obtained on the total gas-to-dust mass 
ratio jMittal et al.|2011| >. 

While the SPIRE and PACS 100 /mi observations were made 
close in time to the IRAM sub-mm observations (2010-08-24 and 
2010-09-09), the PACS 70 /im observations were made half a year 
later (2010-03-14), and so the underlying contribution from the 
AGN may have varied because of the variability in the AGN out- 
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Figure 9. Dust emission. Left: PACS 100 /im image at a resolution of 7" showing the extended dust features. Right: PACS images convolved with a gaussian 
to yield a FWHM of 12". The 70 jjm image is shown in blue, the 100 /im image in green and the 160 /im image in red. 



Table 5. Dust and AGN contributions based on the best-fitting model given 
in Tableg] 



Wavelength /im 


Dust (mJy) 


AGN (mJy) 


20 


2349 ± 221 


88 ± 62 


25 


3384 ± 163 


112 ±74 


30 


3890 ± 282 


137 ± 85 


60 


6772 ± 295 


289 ± 143 


70 


7840 ± 389 


341 ± 160 


100 


8513 ± 775 


502 ± 202 


160 


5532 ± 746 


833 ± 265 


250 


2493 ± 369 


1349 ±323 


350 


1175 ± 171 


1939 ± 354 


500 


485 ± 67 


2850 ± 363 


1153 


50 ±6 


7023 ±714 


1303 


35 ±4 


7389 ± 398 


1428 


27 ± 3 


7673 ± 400 


2710 


4.2 ± 0.5 


10008 ± 396 


2914 


3.4 ± 0.4 


10313 ± 394 


3019 


3.1 ±0.3 


10466 ± 392 


20675 


0.0 ± 0.0 


23231 ± 130 




10.5 11.0 11.5 12.0 12.5 13.0 

Log(^) (Hz) 



O PACS (2010) 


A SPIRE (2010) 


IRAM (2010) 


X UMRAO(2010) 


+ PLANCK (2009) 


□ L'MRAO (2009) 


% SCUBA (1998) 


% WMAP(2002) 


SPITZRR 


# IRAS 


+ FIT 


— — Synchrotron 


— — Cold Dust 


— — Wurm Dusl 


Synchrotron+Dust 



put. According to the best-fit model, this is not an issue of concern 
since the AGN contribution to the total flux at 70 jim is small. This 
was also verified by discarding the 70 yum data point and refitting 
the data. 



4 WHAT ARE THE PREVAILING HEATING 

MECHANISMS IN THE INNER 4 KPC OF NGC 1275 

The filamentary nebula extends out to ~ 50 kpc from the core. The 
excitation mechanisms in the outskirts of the galaxy are very likely 
different from those prevailing in the core. This is evident also from 
the reversal of the relative strengths of [Cn] and [Oi] line emission 
(section[3j2j with cluster-centric distance. The core of NGC 1275, 
where we expect photoionization from stars and AGN to play an 
important role, needs to be modelled separately from the filaments. 
To this end, we conducted simulations using the radiative transfer 
code CLOUDY perland et al.|1998| >, the main goal of which was to 



Figure 10. The spectral energy distribution for NGC 1275. The dust emis- 
sion dominates at FIR wavelengths, whereas the AGN synchrotron emis- 
sion dominates at sub-mm wavelengths. The points overlaid with black di- 
amonds were used to find the best-fit parameters of the model. 



determine whether or not an additional form of heating is required 
to reproduce some of the observed emission lines emerging from 
the core. For example, in the case of NGC 4696, the BCG of the 
Centaurus cluster, the [Cn]/L F1R and Ha/fCn] ratios clearly call for 
another heating component, in addition to photoionization l jMittal| 
|etal.|2011| >. 

From the point of view of energetics of a cooling plasma, ear- 
lier works, such as Johnsto ne et al.|fl9 87 1; Heck man et al.|fT989| l; 
|Voit et al.| ( fT994t ; |Jaffe"& Bremer| p997) ; |Donahue et al.|p000l >, 
have shown that the observed filamentary emissions are far too lu- 
minous to be just due to the recombination phase of gas cooling 
from the ICM. The observed luminosities imply far too much cool- 
ing, with mass deposition rates inconsistent by orders of magnitude 
with the star formation rates. While such calculations have been 
performed for emissions such as Ha and molecular hydrogen lines, 
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Figure 11. The total far-infrared luminosity (8 /im to 1000 /im) vs. the 
100 nm luminosity for seven of the 11 BCGs in the Herschel cool-core 
BCG sample. NGC 1275 is denoted as a black cross. 



none of the current cooling calculations simulate the FIR emis- 
sion lines. This is because the existing cooling calculations (that 
use codes like cloudy, apec, mkcflow) usually stop when the gas 
reaches a temperature of ~ 10 4 K, not low enough to produce the 
FIR lines under consideration. Hence it is presently not possible to 
estimate the luminosity of the FIR lines from a simple cooling-flow 
plasma. 

We modelled a composite cloud comprising a photo- 
dissociation region (PDR) adjacent to an ionized (Hn) region. For 
the equation of state we assumed a constant gas pressure through- 
out the cloud. As input, we provided the normalization for an old 
stellar population (OSP), 77osp, the normalization for a young stellar 
population (YSP), GO (in Habing units), the starting hydrogen den- 
sity, n (in cm -3 ), and the hydrogen column density, Nh (in cm -2 ). 
For further details on the role played by these input parameters, we 
refer the reader to Mittal et al. (201 1). The parameters, along with 
the range explored, are given in Table [8] Based on some initial in- 
vestigations carried out to assess the effect of an OSP component 
on the output emissions, we concluded that the [Cn]/L F i R ratio can 
be better reconciled if an OSP component is included. An old gen- 
eration of stars produces optical and UV photons that contribute to 
the heating of dust but not so much to the photoelectric heating of 
the gas. In the absence of an OSP component, the [Cii]/Lfir for a 
given GO is higher. The results, however, are not very sensitive to 
the presence of an OSP, except for low GO (< 100 Habing) but a 
low GO is not consistent with other observed ratios, as we will show 
below. 

We used the starburst99 stellar synthesis library ( |Leitherer| 
|et al.|1999| > to model both the OSP and YSP. The main input pa- 
rameters to the starburst99 simulations are the stellar mass and age 
(continuous or starburst), and the initial stellar mass function. We 
did not explore the range of possible values these input parameters 
can take - that is a subject of another detailed ongoing study (Mit- 
tal et al. 2012, in prep.). For the purpose of conducting the cloudy 
simulations described here, we fixed the age and mass of the OSP 
to ~ 10 Gyr and 5 x 10" M G , and the age and the mass of the YSP 
to ~ 2 Myr and 5x10' M Q . The SED corresponding to a starburst 
2 Myr old is indistinguishable from the one corresponding to a con- 



Table 6. The default ISM gas phase chemical composition used in cloudy 
simulations. The abundances are given relative to H. 



Element 


log| abundances (ISM) 


He 


-1.0088 


C 


-3.8222 


N 


-4.2010 


O 


-3.7171 


Ne 


-4.1319 


Mg 


-5.1215 


Si 


-5.7222 


S 


-4.7113 


CI 


-7.2218 


Ar 


-5.7716 


Fe 


-6.4218 



tinuous star formation scenario, where the oldest stars are assumed 
to be 2 Myr. This is because the emission in both the cases is dom- 
inated by young stars that are about 2 Myr old. Hence the young 
population of stars need not have formed at the same time. 

The elemental abundances were initially set to their default 
ISM values given in Table [6] However, X-ray observations made 
with Chandra indicate a slight drop in the central metallicity with 
the average value around 0.6 Z Q ( |Sanders & Fabian 2007). With 
this in mind, we also conducted simulations with a lower metal- 
licity, which, as shown below, indeed fit the observed ratios better. 
We did not include any polycyclic aromatic hydrocarbon (PAH) 
grains since the mid-infrared spectra do not contain any PAH fea- 
tures ( |Weedman et al.|2 005 ). Although several BCGs show strong 
PAH features in their IR spectra indicative of star formation ( |Don-| 
|ahue et al.|20TT] >, there are a few such as NGC 1275 and NGC 4696 
( |Kaneda eTa l. 2005 ) that do not. PAH molecules are small in size 
and can easily be destroyed through physical sputtering or thermal 
evaporation jDwek & Arendt|199 2 ; Micelo tta etal.|2011J . 

The absolute and relative strengths of the FIR and Ho- lines 
are very important diagnostics of the various heating contributors 
in cool-core BCGs. These emission lines and their ratios relative 
to [Cn] are listed in Table [7] We note that there is a slight off- 
set of ~ 3" between the peak of the [Cn] emission and the radio 
core emission (the latter coincides with the peak of Ho- emission 
to within ~ 0.5")- This offset is on the order of the lcr pointing 
inaccuracy of the PACS spectrometer and unlikely real. However, 
we estimated the FIR line fluxes for both the cases - (a) assuming 
the offset is not real and (b) assuming the offset is real. These cases 
are referred to as "[Cn] core" and "radio core", respectively. The 
two cases have only marginal differences in their FIR line fluxes 
and are compared only for the purpose of illustrating the level of 
uncertainty in the estimated line parameters. In the following, we 
used the "[Cn] core" as the nominal case and used the other set of 
line fluxes to derive uncertainties on the line ratios. 

The total FIR flux associated with the core can not be directly 
estimated from fitting the SED due to insufficient resolution avail- 
able beyond ~ 160 /jm. We instead used the 100 flux as a proxy 
for the total luminosity. Shown in Fig.[TT]is the total FIR luminosity 
in the range 8 /im to 1000 /zm vs the 100 /jm luminosity for seven 
of the 1 1 BCGs in the Herschel cool-core BCG sample, including 
NGC 1275. We considered only those here because their FIR lu- 
minosity can be determined with least uncertainty. For NGC 1275 
we set Lfir equal to the total dust luminosity. This relation is be- 
ing investigated for the whole sample separately (Oonk et al. 2012, 
in prep.). There is a clear correlation between the total FIR lumi- 
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Table 7. A comparison of FIR and optical emission lines in the core region. Given in the last column are the ratios with respect to [Cn]. The Ho- flux has been 
corrected for the Galactic extinction. We also give the Ha flux after correcting for the internal extinction assuming E(B-V)=0.37 in brackets. 



Region RA 


Dec 


Aperture Radius 


Line 


Velocity 


Flux 


Ratio 






(arc sec) 




km s — ^ 


(1U erg s cm ) 




Radio Core 03hl9m48.16s 


+41d30m42.1s 


1 1 


[Cn] 


11 ±3 


665.1 ± 6.3 


1.00 








[Oi] 


27 ± 4 


1310.5 ± 18.4 


1.97 








[Oib] 


50 ±9 


81.2 ±2.1 


0.12 








[Nil] 


-2 ±7 


53.7 ± 1.1 


0.08 








[Om] 


-26 ± 15 


59.0 ±3.1 


0.09 








Ha 




1810(4185) 


2.7 (6.3) 








LmR 




~ 5.9 x 10 5 


~ 800 


CII Core 03hl9m48.01s 


+41d30m44.9s 


11 


[Cn] 


8± 3 


724.8 ± 7.3 


1.00 








[Oi] 


23 ±4 


1320.2 ± 17.3 


1.82 








[Oib] 


47 ±8 


84.3 ± 2.0 


0.12 








[Nil] 


-4± 8 


58.1 ± 1.2 


0.08 








[Om] 


-30 ± 14 


66.1 ± 3.0 


0.09 








Ha 




1810(4185) 


2.5 (5.8) 








LmR 




~ 5.9 x 10 5 


~ 700 



nosity and the 100 /im luminosity, such that L FIR oc L\^ m . The 
trend is also observed in the flux-flux plane (not shown) and so the 
correlation is not spuriously induced due to the common depen- 
dence of the two quantities on redshift. Using this correlation and 
the 100 fim core flux, the total FIR flux for the core may be es- 
timated. The contribution from the AGN synchrotron emission at 
100 j-im is negligible (< 6%, see Table[5} and may be ignored. 

The measured 100 yum flux within the core aperture of 11" 
radius is (5608±7) mjy. Even though some of the flux may fall 
beyond the chosen aperture, necessitating an aperture correction, 
the factor is very small for the given radius. We determined the re- 
quired aperture correction factor by assuming that the FIR emission 
correlates with the Ho- emission. Assuming the WIYN Ho- image 
to represent the true model for the FIR surface-brightness distribu- 
tion, we convolved the Ho- image with the PACS 100 /im PSF and 
estimated the aperture correction from the ratio of the flux from the 
smoothed to the unsmoothed image. The ratio is close to unity and 
therefore we did not apply any aperture correction. The measured 
100 yum flux yields a FIR (8 fim to 1000 fim) flux that implies a 
£fir/[Cii] ratio ~ 700 to 800. 



4.1 Best-fit energy model for the core 

The Herschel FIR coolants serve as strong constraints for evaluat- 
ing the physical parameters of the ISM in the BCG. Fig. |12| shows 
[Oi]/[Cii] versus n for different values of the photon parameter, GO. 
This plot can be understood in terms of the critical densities for 
[Cn] and [Oi]. This is the density at which the probabilities of col- 
lisional de-excitation and radiative de-excitation are equal ( |Oster-| 
|brock &"Fe rland 20061. Above the critical density the level pop- 
ulations may be considered to be in local thermal equilibrium as 
the level populations become dominated by collisions. The criti- 
cal density is lowered by a factor that is roughly the optical depth 
of the line due to photon trapping when the line is optically thick. 
Fig. [12] shows that, for n < n cr , the intensity of an optically thin 
line increases linearly with n and for n > n a the intensity becomes 
independent on n. This result is fairly insensitive to details such as 
the cloud energy source, as Fig. [12] shows. While n CI for [Cn] is 
~ 3 x 10 3 crrT 3 , it is substantially higher for [Oi] ~ 5 x 10 5 cnT 3 . 
Therefore, at a fixed column density and photon parameter, GO, for 
high densities 3 x 10 3 citT 3 < n < 5 x 10 5 cnT 3 the ratio, [Oi]/[Cii], 



continues to increase with density. This makes [Oi]/[Cn] a sensitive 
probe of density if the lines are optically thin. The individual trends 
displayed by the [Oi] and [Cn] lines in Fig.[l2]are complicated by 
the fact that for a given column density these lines become optically 
thick for certain combinations of n and GO. 

Similarly, Fig.[T3]shows [Nn] and [Om] versus n for different 
GO and for two different column densities. [Cn] is produced in both 
ionized (for e.g. Hn regions) and neutral media (for e.g. PDRs), 
however, [Nn] and [Om] are produced only in ionized region due 
to their higher ionization potentials. This has the effect that [Nn] 
and [Om] are produced profusely at the surface of the cloud facing 
the ionization source and reduce with the depth into the cloud. This 
explains the drop in both [Nn] and [Om] intensities relative to [Cn] 
with Nh- This makes [Nii]/[Cii] and [Om]/[Cn] a sensitive probe of 
the cloud depth or hydrogen column density, N H . 

Of the Nh range investigated, the optimal cloud depth re- 
quired to reproduce the observed ratios is 10 23 cm 4 . Unless GO 
is high (> 10000 Habing), a higher N H does not affect the pre- 
dicted ratios, implying that the cloud becomes radiation bounded. 
In Fig. [14] we present the modelled ratios as function of n and GO 
for N H = 10 23 cm" 2 , assuming stellar photoionization only. The dif- 
ferent sets of curves correspond to the observed lower- and upper- 
limits of [Cn]/L F1R (solid red), [Oi]/[Cii] (dotted blue), Ha/[Cn] 
(dashed green), [Nii]/[Cii] (dash-dotted orange), [Oib]/[Cii] (dash- 
dot-dotted pink) and [Om]/[Cn] (long-dashed gray). Note that the 
observed Ho'/[Cn] ratio has a wide range due to the uncertainties 
mentioned in section [2~2l The FIR ratios, in contrast, are much bet- 
ter constrained. Consequently, our strategy was to first try and re- 
produce the FIR line and continuum ratios and then use Ha/[Cn] to 
check for consistency against the best-fit model. 

The left panel of Fig.[l4]assumes ISM abundances, as adopted 
in |Ferland et aL] < f2008l l. While [Oi]/[Cii], [Om]/[Cn] and [Om]/[Cn] 
converge at n ~ 1000 cirT 3 and GO ~ 1000 Habing, the predicted 
[Nn]/[Cn] for this combination of n and GO is less than the observed 
ratio. The middle panel of Fig. [14] assumes the nitrogen abundance 
is twice the ISM value, yielding a much better convergence between 
the FIR ratios. Curiously, the best-fitting («,G0) for NGC 4696 in 
the Centaurus galaxy cluster can also be reconciled with the ob- 
served [Nii]/[Cii] only if the nitrogen abundance is increased by 
a factor of two over the assumed overall metallicity (Mittal et al. 
|201 1| >. In the case of NGC 4696, the need for an overabundance 
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Figure 12. [Oi]/[Cn] sensitivity to the gas density. The simulated [Cn] intensity in erg s" 1 cm" 2 (upper left panel), [Oi] intensity in erg s" 1 cm" 2 (upper 
right panel) and [Oi] to [Cn] ratio (lower panel) as a function of gas density, n, and the normalization of the young stellar population (the photon parameter), 
GO. The column density has been fixed to Nyj= 10 23 cm" 2 . The two vertical dashed lines represent the critical densities, n a - = 3 X 10 3 cnT 3 for [Cn] and 
n a = 3 X 10 3 cnT 3 for [Oi]. The horizontal solid line corresponds to the observed ratio of [Oi] to [Cn] for the inner 4 kpc region. 



Table 8. Input parameters for the cloudy simulations. 



Parameter 


Symbol 


Input Range 


Likely Values 


Total Hydrogen Density (cm -3 ) 


n 


10 to 10 6 


500 to 700 


FUV Intensity Field (Habing*) 


GO 


1 to 10 6 


700 to 900 


Hydrogen Column Density (cnT 2 ) 


N H 


10 19 to 10 26 


> 10 23 


Metal licity 


Z 




0.6 


Nitrogen abundance (relative to Z) 


Ze(N) 




2 


Normalization for the OSP* (10" 16 erg s" 1 cm" 2 Hz" 1 ) 


7?OSP 




52.54 



a 1 Habing = 1.6 X 10" 



t -1 -2 

erg s cm , 



OSP: Old Stellar Population 



in nitrogen is confirmed by independent observations in the opti- 
cal and X-ray bands. An overabundance of nitrogen could imply an 
abundance of intermediate-mass asymptotic giant branch (AGB) 
stars with pre-main sequence mass between 3 M Q and 8 M B , which 
may convert dredged-up carbon in their outer shells into nitrogen 
(Wood et al.|1983l|Iben|1975l >. We are conducting such detailed ra- 
diative transfer modelling of all the cool-core BCGs in the Herschel 
sample and it will be interesting to see if this is a generic feature of 
galaxies at the center of cooling-flows. 

As mentioned in section]?] X-ray observations suggest a lower 
metallicity in the core of Perseus (< 20" ~ 7 kpc). Hence, in the 
right panel of Fig.[l4] we show the predicted ratios for a metallicity 
equal to 0.6 times the ISM value. The FIR ratios coincide nicely 
at n ~ 650 cm" 3 and GO ~ 800 Habing. We can now derive con- 
straints on the physical scale of the cloud by comparing the [Cn] 
flux emerging from the modelled cloud and the observed [Cn] flux. 
For the range of n and GO given in Table[8] the emergent [Cn] flux 



cm 2 to 3.7 x 10 3 erg s 1 cm 2 . 



is in the range 3.2 X 10 erg s" 
The size of the cloud that produces [Cn] flux at the distance of Earth 
between 665.1 x 10" 15 erg s" 1 cm" 2 and 724.8 x 10" 15 erg s" 1 cm" 2 
is > 1 kpc, yielding a volume filling factor of 1.7 x 10" 2 . 

The main conclusion from modelling the ISM associated with 
the core of NGC 1275 is that a stellar-like ionizing source with 
GO ~ 1000 Habing is capable of reproducing the observed FIR 
emission and Ha emission. No additional form of heating is re- 
quired to explain these observations. 



4.2 Consistency with optical flux ratios 

According to previous studies iJohnstone & Fabian 1988, |Don-| 



|ahue & VoitlfTWT] [Johnstone et al.||2007| |Ferland et al^OO"^ 
pure stellar or a pure active nucleus origin as the source of the 
optical filaments has difficulty explaining certain optical flux ra- 
tios. This may be true for the filaments but we suspect that stel- 
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Figure 13. [Nii]/[Cii] (upper row) and [Om]/[Cii] (lower row) versus density, n. Shown are the ratios for two different column densities, Nh= 10 19 cm -2 (left) 
and Nh= 10 23 cm -2 (right). [Nn] and [Om] are produced in abundance at the surface of the cloud facing the ionization source and deplete with the depth of 
the cloud and, hence, sensitive probes of Nh_. The different points have the same meaning as in Fig. |12| The horizontal lines represent the observed ratios in 
the inner 4 kpc region. 
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Figure 14. Simulations based on a photoionization model containing only stars. Shown are contours corresponding to the lower and upper observed limits 
(given in square brackets below) for L mR /[Ca] (solid red) -» [700:800]; [Oi]/[Cii] (dotted blue) -> [1.79:2.00]; Ha/[Cn] (dashed green) -> [2.5:6.3]; [Nii]/[Cii] 
(dashed-dotted orange) -» [0.078:0.082]; [Om]/[Cn] (dash-dot-dotted pink) -» [0.113:0.125] and [Om]/[Cn] (long-dashed gray) -> [0.083:0.095]. For this 
simple model, the regions allowed by the data, between the paired lines, would be expected to overlap at a single Habing flux (GO) and cloud density («). Left: 
The element abundances are fixed to their default ISM values. Middle: Nitrogen overabundance is fixed to two. Right: The element abundances are fixed to 0.6 
times the ISM value and nitrogen overabundance to two. The best-measured ratio, Lfir/[Cii], overlaps for this model. 



lar photoionization becomes increasingly important with decreas- 
ing distance to the core of the galaxy. Although fitting line and 
continuum ratios measured in different bands across the electro- 
magnetic spectrum (submm, IR, optical, UV) simultaneously is be- 
yond the scope of this paper, we compared a subset of the opti- 
cal ratios measured by Johnstone & Fabian ( 1988} to our best-fit 



model predictions. This comparison is shown in Table [9] We give 
a range for the observed ratios since the ratios display a strong 
radial gradient. 'Model F refers to the best-fit energy model de- 
scribed above. While the predicted values of the [Nn]/i6583/Ho' 
and [Sn]/16731/Hcj' ratios are consistent with the observed val- 
ues, [On]/i/!3727, 3729/Hy3 is lower by a factor of a few and 
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Table 9. The optical emission line ratios in the central 4 kpc region of 
NGC 1275 (Johnstone & Fa bian|1988) . The columns are (1) the species, 
(2) the measured ratios (3) the predicted ratios based on the best-fit model 
(4) the predicted ratios based on a slightly modified model. The observed 
ratios are from the observations off to the south-east side of the nucleus to 
avoid the high-velocity system. 



Ratio 



Value 



Model I 



Model II 



[On]/U3727, 3729/H/3 3.5-8 
[Ora]i5007/Hj8 1.2-1.8 

[Oi]/16300/Ha+[Nii]i6583 0.06 - 0.14 
[Nn]i6583/H(»- 0.8-1.4 
[Sny6731/Ha 0.2-0.7 



1-2 1-2 

0.1-0.2 0.7-2 

0.004-0.007 0.008-0.01 

0.6-0.9 0.6-0.9 

0.6-0.8 0.6-0.8 



[Oiii]^5007/HvS and [Oi]/l6300/[Nii]^6583+Ha are lower by about 
an order of magnitude than the observed values. We devised another 
model, referred to as 'Model II', obtained by slightly modifying the 
best-fit model by decreasing the age of the input young stellar pop- 
ulation by an order of magnitude. Using Model II we are able to re- 
produce some of the optical ratios, specifically, the ratio [Om]/H/?. 
There is also an increase in the ratio, [Oi]/l6300/[Nii]/i6583+Ha', 
by a factor of two. The modified model does not impact the best-fit 
model parameters constrained by the relative strengths of FIR and 
Ho- lines (Table [7} . This is due to the fact that lowering the age of 
the oldest stars in the YSP alters the SED only below 1000 A, such 
that [Om]/I5007 and [Oi]/16300 are affected. FIR lines, on the other 
hand, are mainly produced by gas heated by lower-energy UV and 
optical photons. 

From this exercise we conclude that photoionization from 
stars is the predominant source of energy in the cold gas clouds in 
the core of NGC 1275. It may be that in order to obtain consistency 
with the complete set of observed flux ratios at different wavebands, 
a composite model that includes a stellar component and one or 
more additional forms of heating is required. Particle-heating mod- 
els are usually very hard to distinguish from slow shocks (e.g. Fer- 
|land et al.|2008||Farage et al.|2010||Rich et al.|201l) . Even though 
the work of |Ferland et al.| ( |2009[ l showed that ionizing particles 
are the likely heating agents in the filaments of NGC 1275, the 
particle heating model was only marginally more successful than 
heating by dissipative MHD wave energy. The core of NGC 1275 
displays large-scale motions, which may very well entail shocks 
(although strong shocks have been ruled out by X-ray data, Fabian 
et al. 2003, Sanders & Fabian 20071. Similarly, due to the proxim- 
ity to the AGN, particle heating (cosmic rays) may be significant 
as well. Hence additional heating required to explain the FIR and 
optical line ratios in entirety may either manifest itself in the form 
of weak shocks or particle-heating or both. 



4.3 Star formation 

Far-ultraviolet (FUV) emission is a direct indicator of mor- 
phological and spatial extent of recent star-formation sites. 
In Fig. [15] we show the FUV emission associated with the 
BCG of Perseus. This image was created using two datasets 
from the HST archive - (1) Space Telescope Imaging Spectro- 
graph (STIS) FUV-MAMA/F25SRF2 data (proposal ID: 8107, PI: 
C. O'Dea) and (2) Advance Camera of Surveys Solar Blind Chan- 
nel (ACS/SBC) FU V-M AM A/F 1 40LP data (proposal ID: 11207, 
PI: R. O'Connell). The first dataset consisted of a single pointing 
only (centered on the core). The second dataset, however, consisted 
of 8 different pointings, designed to map the bright UV filaments 



extending NW and SE from the core. We obtained the single flat- 
fielded, dark-subtracted exposures from the HST archive for both 
the datasets and combined them using the iraf task, multidrizzle. 

The filaments NW of the core have the same spiral morphol- 
ogy as the HVS. This spatial correlation seems to suggest that the 
star clusters in the filaments are associated with the foreground 
galaxy (e.g. Keel & White 2001). There is also evidence to the 
contrary, suggesting that the star clusters belong to the LVS (e.g. 
|Goudfrooij|l995||Brodie et al.|1998| >. 

For the present study, we will focus only on the FUV emis- 
sion originating from the core region (see Fig. [16} , which, due 
to proximity and symmetry, is expected to be associated with the 
LVS. The integrated flux-density within the core is (5.3 ± 0.1) x 
10~ 15 erg s cnT 2 A -1 at 1456 A; however the FUV emission 
close to the center is likely to be associated with the AGN. For 
the purpose of quantifying FUV emission from stars, we deter- 
mined the FUV emission in the annulus (0.5 kpc - 3.8 kpc) as 
(3.6 ± 0.1) x 10~ 15 erg s -1 cirT 2 A , and correcting for the Galactic 
extinction results in (1.2 ± 0.3) x 10~ 14 erg s -1 cirT 2 A -1 . There 
is a further correction factor that needs to be applied in order to 
account for the dust internal to the ISM of NGC 1275 (see sec- 
tion |2.3} . Assuming an internal reddening of E(B - V) = 0.37, we 
obtain a FUV flux-density of (1.9 ± 0.5) x 10~ 13 erg s -1 cm" 2 A" 1 . 
This can be directly compared to the flux density expected from a 
synthetic spectrum of a young stellar population, such as the one 
used as input for the cloudy simulations. The expected flux density 
was determined by convolving the redshifted synthetic spectrum 
with the bandpass of the FUVMAM A/F25 SRF2 filter on STIS us- 
ing the iraf tool synphot. A synthetic spectrum corresponding to 
an instantaneous starburst containing 5xl0 9 M Q and 2 Myr in age 
predicts a flux density of 1.7 x 10~" erg s -1 cirT 2 A -1 . Holding 
the age constant and scaling down the mass so that the flux density 
agrees with the extinction-corrected observed flux density implies 
a SFR of 27 M e yr 1 . This is similar to the SFR ~ 24 M Q yr 1 de- 
rived in section [33] (Table |4} from the FIR measurements. This is 
also in very good agreement with the SFR ~ 25 M Q yr~' derived 
by |Norgaard-Nielsen et al.|(T99(3) from the IUE (International Ul- 
traviolet Explorer) spectroscopic data. 

Since the star formation rate traced via FIR data is expected to 
be most accurate (due to negligible extinction at FIR wavelengths), 
and the FUV rate matches it, the average internal reddening over 
the inner 4 kpc radius region is probably close to that suggested 
by the observed Galactic-extinction corrected Balmer decrement 
(4.07) in a region 18" SW of the nucleus (section [2~3| >. On the other 
hand, the SFR derived above using FUV data is for the core re- 
gion, whereas the SFR derived in Table [4] is for the whole galaxy. 
For example, Canning et al. (20101 obtained an optical star forma- 
tion rate of about 20 M 3 yr~' over the Blue Loop region, which is 
outside the core region. Although there is hardly any FIR emission 
detected in the Blue Loop region, it is possible that it contributes 
to the total SFR derived from FIR measurements to some extent. 
This could imply that the average internal reddening over the core 
is lower than assumed, making the FUV-derived star formation rate 
lower as well. Hence E(B — V) ~ 0.37 is an upper-limit. Note, 
the above calculation is based on two assumptions: (1) the dust in 
NGC 1275 follows the same extinction law as the Milky Way and 
(2) the Balmer line ratios are close to the case-B recombination 
values. 

It is interesting that if we use the dust attenuation derived from 
the Balmer decrement and apply it at smaller (UV) wavelengths, the 
FIR- and UV-derived SFRs are consistent with each other. This is in 
contrast to studies that have indicated that the reddening values di- 
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Figure 15. Combined HST SBC/F140LP and STIS/F25SRF2 FUV image 
of NGC 1275. The red circle indicates the central 4 kpc core region. 



verge between nebular line emission and UV continuum by a factor 
of two (e.g. |Calzetti|1997|[Buat et al.|2002) , such that nebular lines, 
like Ha and H/3, are more attenuated than the stellar continuum 
emission. This discrepancy is usually attributed to an uneven dis- 
tribution of dust in front of stars and ionized gas, with the latter be- 
ing more closely associated with the dust than the former. Calzetti 
( 1997) suggests using the relation E(B- V) sb „ = 0.44 x E(B- V) gas . 
However, |Garn et aT| (20 1 ) argue that the while the continuum at 
longer wavelengths (optical) may be probing older stellar popula- 
tion residing in less dusty environments, and hence may not be co- 
spatial with the ionized gas, the continuum at shorter wavelengths 
(UV) should be more closely tied to the young stars producing ion- 
ized gas. Hence the correction factor (0.44) for FUV emission may 
be higher than suggested by the above relation. 

Using the relation suggested by Calzetti| {1997[ >, |Buat et aT| 
(20021 derived ^uvy=2000) = 1.6 X Aa a . Instead of applying the at- 
tenuation at FUV wavelengths by extrapolating the reddening value 
calculated from the Balmer decrement, if we apply the above cal- 
ibration we obtain a SFR of 6 M Q yr~'. This is a factor of four 
lower than the FIR-derived SFR. This implies that in the case of 
NGC 1275, the dust covering factor for the ionized gas resulting 
in the Balmer lines and the FUV emission must be rather homoge- 
neous, such that the attenuation tracked by the ionized gas is closely 
related to that of the FUV emission. 



5 DISCUSSION OF THE OUTER FILAMENTS 



Ferland et al. ( 2009, F09 hereafter) conducted a thorough analysis 
of the Horseshoe region in NGC 1275. They used the infrared and 
optical line intensities to distinguish between two types of heating: 
extra heating, as would be produced by passing shocks, and heating 
by energetic particles such as cosmic rays. Extra heating refers to 
heating by dissipative MHD waves, and it is assumed that this kind 




Figure 16. A zoom-in FUV image of the center of the galaxy, the red circle 
indicates the central 4 kpc core region. 



of heating increases the thermal energy of the gas only. Hence the 
local gas kinetic temperature dictates the collisional processes that 
are energetically possible. Energetic particle heating refers essen- 
tially to ionizing particles, and, on the other hand, is not only capa- 
ble of depositing thermal energy via elastic collisions with the par- 
ticles of the ISM but can also ionize a neutral medium via creation 
of a population of suprathermal secondary electrons. The physics 
of molecular gas exposed to such heating sources is described by 
|Ferland|pOTTat and |Ferfand|(201 lb| >. 

The model in F09 assumes that the gas in the filaments is 
in pressure balance with the surrounding X-ray-emitting intraclus- 
ter gas. Based on the X-ray measurement of the electron pres- 
sure of 0.128 keV cirr 3 of the hot gas surrounding the horseshoe 
by Sanders & Fabian (2007), F09 assumed a constant pressure of 
nT = 10 65 cnT 3 K. The simulations conducted in F09 are fur- 
ther based on the assumption that the filaments are composed of 
cloudlets with a range of density, n, and temperature, T, but which 
have this single pressure, nT. This assumption is motivated by the 
fact ionic, atomic, and molecular emissions are all observed in the 
filaments. This implies that different phases of gas occupy a tele- 
scope beam, even at HST resolution. In other words, there is ob- 
servational evidence that both dense molecular and diffuse ion- 
ized emission arise from spatially coincident regions. Because of 
the constant pressure assumption, low density cloudlets have high 
temperature and produce emission from ionized gas, while dense 
clouds are cold and account for the molecular component. 

The formalism adopted in F09 uses a cumulative filling factor 
f(n), which is a powerlaw in density, as the weighting function to 
co-add clouds of different densities. This factor describes the frac- 
tional volume filled with gas with density n or lower. The spectra 
for various emission lines, in particular the infrared Hi and opti- 
cal H i emission lines, are determined using cloudy using a range 
of electron densities, temperatures, and non-radiative heating rates. 
The emission for a given line is integrated over the ensemble of 
clouds and then compared to observations. 

|Ferlan d et al. ( 2009) found that both the forms of non-radiative 
heating, extra heating and energetic particle heating, match the op- 
tical and infrared observations to within a factor of two for the ma- 
jority of the lines. There are a few discriminant lines, such as the 
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Table 10. A comparison of [Cn], [Oi] and Ha line fluxes in the Horseshoe knot, the SW1/SW2 knots and the Blue Loop knots. Given in the last column are 
the ratios with respect to [Cii]. [Oi] is not detected in any of the extended filaments - we give the 3<r upper-limits. 



Region 


RA 


Dec 


Aperture Radius 
(arcsec) 


Line 


Velocity 
km s 


Flux 

(10~ 15 erg s _1 cm- 2 ) 


Ratio 


Horseshoe knot 


03hl9m45.15s 


+41d31m33.1s 


11 


[Cn] 


79.5 


26.78 ± 0.26 


1.00 










[Oi] 




<44 


< 1.64 










Ho 




36 (88) 


1.34 (3.27) 


Southwest knots 


03hl9m45.48s 


+41d29m49.4s 


11 


[Cii] 


3.3 


40.67 ± 3.55 


1.00 


(Southern Filament) 








[Oil 




< 34 


< 0.84 










Ho 




87 (211) 


2.14(5.19) 


Blue Loop knots 


03hl9m50.72s 


+41d30m07.2s 


11 


[Cn] 


-224.5 


35.19 ± 1.87 


1.00 


(Southeast Filament) 








[Oil 




< 30 


< 0.85 










Ho- 




73 (177) 


2.07 (5.03) 



optical emission lines He i A 5876 A, [Ne in] A 3869 A, and the 
infrared emission line [Ne n] A 12.81 /jm, which show a few or- 
ders of magnitude difference and indicate that ionizing particles are 
responsible for heating and ionizing the gas. |Fabian et al.| |20TTJ 
argue that the surrounding hot intracluster medium is the source 
of the ionizing particles, rather than true cosmic rays as would be 
found in the adjoining radio lobes. 

All of the predicted lines used to compare with observations 
were optically thin, so their intensity relative to similar forbidden 
lines has no dependence on cloud column density. The F09 model 
also made predictions for the far-infrared emission lines, such as 
those observed by Herschel. In Table |10| we list the detected Her- 
schel emissions in the three regions of the extended filaments - 
the Horseshoe region, the southwest (SW) knots. We also give the 
3-<x upper-limits for the non-detections. These lines may be used 
to distinguish between the two heating scenarios. In particular, the 
predicted [Oi]/[Cii] ratio is ~ 3 for the extra heating and ~ 21 for 
the energetic particle heating. The observed [Oi]/[Cii] ratios, on the 
other hand, have an upper-limit of 1.64 in the Horseshoe knot and 
0.85 in both the SW and Blue Loop knots. This discrepancy was 
first pointed out by Mittal et al. (201 1 1 in a different context. 

There are two plausible explanations for why the F09 model 
fails to reproduce the Herschel observations. The first reason is tied 
to the critical densities, « cr , of [Cn] and [Oi] gas (see Fig. [72}. Since 
the F09 model assumes a constant pressure of nT = 10 6 5 cm -3 K, 
the FIR lines which are produced at low temperatures (~ 100 K to 
a few hundred K ) correspond to high densities (3 x 10 4 citT 3 > n > 
5 x 10 3 citT 3 ). Such high densities lead to high [Oi]/[Cn] ratios due 
to the reasons given in section [4~T| The model could be brought into 
agreement with the observations by postulating a large component 
of low density gas, which would strongly emit the Herschel lines 
with the observed ratio. F09 note, in section 7, point 7, that large 
reservoirs of cold gas could be present yet not detected with the 
selection of lines they had available. 

The second likely cause for the incompatibility of the F09 
model with the [Oi]/[Cii] ratio is that the F09 model assumes that 
the emission lines are optically thin. If the gas has high enough col- 
umn density, the Herschel lines become optically thick. This is nor- 
mally the case in galactic PDRs ( Tielens & Hollenbach 1985}. The 
line luminosity is no longer determined simply by the product of the 
line emissivity and the volume of the cloud but rather the geometry 
of the emitting cloud, especially its column density, will affect the 
line intensities. F09 note that the observed surface brightness of the 
Ho- line gives the line-of-sight thickness of the cloud Al ~ 0.3 pc. 
This implies a hydrogen column density of N(H) = 10 225 citT 2 . 
If the emission forms in a single cloud with this column density, 



the Herschel lines would be strongly affected but there would be 
little impact on the NIR and optical lines, other than the effects of 
internal reddening. 

Clearly much remains to be learned about the geometry of 
the filaments at sub-HST resolution scales. Both scenarios outlined 
here are consistent with what is known from available observations. 



6 CONCLUSIONS 

We have presented far-infrared Herschel observations of the cen- 
ter of the Perseus galaxy cluster. The brightest cluster galaxy, 
NGC 1275, is surrounded by filaments, previously imaged exten- 
sively in Ha and CO emission. In this work, we report the presence 
of coolants, such as [Cn], [Oi], [Nn], [Oib] and [Om], in addition 
to Hor and CO. All the Herschel lines except [Om] are spatially ex- 
tended, with [Cn] extending up to 25 kpc from the nucleus of the 
galaxy, and cospatial with Ho- and CO. Furthermore, [Cn] shows a 
similar velocity distribution to CO, and the latter has been shown 
in previous studies to display a close association with the Her kine- 
matics. The spatial and kinematical correlation among [Cn], Ha 
and CO gives us confidence that the different components of the 
gas may be modelled with a common heating model. 

The velocity structure inferred from [Cn] observations reveals 
blueshifted components on either side of the nucleus and a ridge of 
redshifted gas passing through the center in the north-south direc- 
tion. A combination of disk rotation on ~ 10 kpc scales and more 
disordered motion of a few distinct clouds at larger radii is one 
plausible scenario consistent with [Cn] (Herschel) and CO data. It 
may also be that outflows from the radio source at the center are 
dredging up cold material away from the line of sight, resulting in 
a ridge of redshifted gas. 

We have detected continuum emission at all PACS and SPIRE 
wavelengths ranging from 70 yum to 500 jim. The Herschel con- 
tinuum data combined with a suite of existing radio, submm and 
IR data, have allowed us to perform a spectral energy distribu- 
tion fitting using a model comprising two modified-blackbody dust 
components, dominant at IR wavelengths, and synchrotron AGN 
emission, dominant at radio and submm wavelengths. The total in- 
frared luminosity (8 yum to 1000 yum) is (1.5 ± 0.05)10" L Q , mak- 
ing NGC 1275 a luminous infrared galaxy (LIRG). The infrared- 
inferred star formation rate is about 24 M Q yr~' . This is comparable 
to the SFR estimated in a core region (~ 4 kpc in radius) using the 
HST-FUV observations, and also XMM-Newton RGS observations, 
which suggest a residual cooling rate of 20 M Q yr~' . 

We have investigated in detail the source of the emissions 
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emerging from a core region 4 kpc in radius. This investigation 
was done by carrying out radiative transfer simulations assuming 
a photoionization model consisting of only old and young stellar 
populations as heating agents. We find that the Herschel and Ha 
emissions can be reproduced by such a model, yielding hydrogen 
density between 500 cirr 3 and 700 cirr 3 and FUV intensity field 
between 700 Habings to 900 Habings. Optical flux ratios indicate 
that a second heating component may be needed; however, stellar 
photoionization seems to be the dominant mechanism. 

We have also detected [Cn] in three previously well-studied 
regions of the filaments: the Horseshoe, Blue Loop and southwest 
knots. The observed upper-limits of the [Oi]/[Cii] ratio are 1 dex 
smaller than predicted by the best-fit model of Ferland et al.| ( |200"9"l >, 
wherein two sources of ionization have been considered individu- 
ally: particle heating and extra-heating (the latter increases the ther- 
mal energy of the gas). The [Cn] line has a excitation potential of 
91 K and a critical density much lower than the lines considered in 
that study. This suggests that the lines are optically thick, as is typ- 
ical of galactic PDRs, and implies that there is a large reservoir of 
cold atomic gas. This has not been included in previous inventories 
of the filament mass and may represent a significant component. It 
seems likely that the model used by Ferland et al. (2009 ) needs to 
be augmented in order to reproduce the observed strengths of the 
various emissions in the filaments of NGC 1275. It can also be that 
a composite model, like the one used for the study of Centaurus 
( |Mittal et al.|201l) , consisting of more than one heating agent, is 
required to achieve compatibility with observations. 
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